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LISTING OF THE CLAIMS 

Prior to calculation of the filing fee, please amend the claims as follows. This 
listing of claims will replace all prior versions, and listings, of all claims as pending in 
this application. 

Clainns 1-10 (Canceled) 

Claim 11 (New): A method of estimating, from data obtained by 
exploration of a zone of a heterogeneous medium, a model representative of a 
distribution, in the zone, of at least one physical quantity, the model being free of a 
presence of correlated noises that may be contained in the data, comprising: 

a) acquiring measurements giving information about physical 
characteristics of the zone by following a predetermined experimental protocol; 

b) specifying a noise modelling operator which associates, with a model 
of each physical quantity, synthetic data that constitute a response of the model, the 
measurements and the synthetic data belonging to a data space; 

c) selecting, for each correlated noise referenced by a subscript] ranging 
from 1 to J, a noise modelling operator which associates a correlated noise with a 
noise-generating function belonging to a predetermined spaced of the noise- 
generating functions (Bj); 

d) specifying a norm or of a semi-norm in the data space; 
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e) specifying a semi-norm in the space of the noise-generating functions 
for which each noise modelling operator establishes substantially an isometric 
relation between the space of noise-generating functions and the data space; 

f) defining a cost function quantifying a difference between the 
measurements on one hand and a superposition of a model response and of the 
correlated noise associated with the noise-generating function on the other hand; 
and 

g) adjustment of the model and of the noise-generating functions by 
minimizing the cost function, by means of an algorithmic method taking advantage of 
isometry of properties of the noise modelling operators. 

Claim 12 (New): A method as claimed in claim 1 1 , wherein: 

a distribution as a function of depth of an acoustic impedance in the medium 
is sought, the correlated noises affecting the data are tube waves each identified by 
parameters characterizing their propagation, the measured data are VSP data 
obtained by means of pickups which detect displacement of particles in the medium 
in response to a localized seismic excitation, the location of the pickups and a 
recording time and time sampling points being defined, and the selected noise 
modelling operator associates the synthetic data with an acoustic impedance 
distribution as a function of an evaluated depth in traveltime and with vertical stress 
measured as a function of time at a depth of the first pickup. 
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Claims 13 (New): A method as claimed in claim 12, wherein: 

the cost function quantifying the difference is a square of the semi-norm of the 
difference in the data space. 

Claim 14 (New): A method as claimed in claim 12, wherein: 

adjustment of the model and of the noise-generating functions is obtained by 
means of a block relaxation method for eliminating unknowns corresponding to each 
correlated noise generating function, the block relaxation method being Implemented 
within iterations of a quasi-Newtonian algorithm for calculation of the model. 

Claim 15 (New): A method as claimed in claim 13, wherein: 

adjustment of the model and of the noise-generating functions is obtained by 
means of a block relaxation method for eliminating the unknowns corresponding to 
each correlated noise generating function, this the block relaxation method being 
implemented within the iterations of a quasi-Newtonian algorithm for calculation of 
the model. 

Claim 16 (New): A method as claimed in claim 12, wherein: 
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numerical calculation of the image of a model by the modelling operator is 
carried out by a numerical solution of a 1 D wave equation for the model, by selecting 
values taken by displacement of the particles at locations of pickups and at 
previously defined time sampling points, and by applying an operator likely to 
compensate for spherical divergence and attenuation effects. 

Claim 17 (New): A method as claimed in claim 13, wherein: 

numerical calculation of the image of a model by the modelling operator is 
carried out by a numerical solution of the a 1D waves equation for the model 
considered, by selecting values taken by the displacement of the particles at the 
locations of pickups and at the previously defined time sampling points, and by 
applying an operator likely to compensate for the spherical divergence and 
attenuation effects. 

Claim 18 (New): A method as claimed in claim 14, wherein: 

numerical calculation of the image of a model by the modelling operator is 
carried out by a numerical solution of the a ID waves equation for the model 
considered, by selecting values taken by the displacement of the particles at the 
locations of pickups and at the previously defined time sampling points, and by 
applying an operator likely to compensate for the spherical divergence and 
attenuation effects. 
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Claim 19 (New): A method as claimed in claim 15, wherein: 

numerical calculation of the image of a model by the modelling operator is 
carried out by a numerical solution of the a 1D waves equation for the model 
considered, by selecting values taken by the displacement of the particles at the 
locations of pickups and at the previously defined time sampling points, and by 
applying an operator likely to compensate for the spherical divergence and 
attenuation effects. 

Claim 20 (New): A method as claimed In claim 12, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing a noise transport equation, and a utilized noise generating function 
has an initial condition along an edge of an observation zone belonging to the space 
(Bj): of support time functions in a give time interval. 

Claim 21 (New): A method as claimed in claim 13, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
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observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 22 (New): A method as claimed in claim 14, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 23 (New): A method as claimed in claim 15, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 24 (New): A method as claimed in claim 16, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
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function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 25 (New): A method as claimed in claim 17, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 26 (New): A method as claimed in claim 18, wherein: 

the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 27 (New): A method as claimed in claim 19, wherein: 
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the noise modelling operator is a finite-difference centered numerical scheme 
for discretizing the noise transport equation, and the utilized noise-generating 
function involved as the function has an initial condition along the an edge of the an 
observation zone belongs belonging to a the space (Bj) consisting of the support time 
functions in a give time interval. 

Claim 28 (New): A method as claimed in claim 12, wherein: 
a semi-norm selected for the data space is: 


1 



and a semi-norm selected for the space of the noise-generating function is: 


1 


\m\B = 




Claim 29 (New): A method as claimed in claim 13, wherein: 


the a semi-nomn selected for the data space is: 
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llwlb = f Ax At + <f] 

\ 1=0 n=0 ^ / 


and the a semi-norm selected for the space of the noise-generating function space 
is: 


Claim 30 (New): A method as claimed in claim 14, wherein: 
the a semi-norm selected for the data space is: 

= (ax At 'y. E («^' + 


and the a semi-norm selected for the space of the noise-generating function space 
is: 
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Claim 31 (New): A method as claimed in claim 15, wherein: 
the a semi-norm selected for the data space is: 

iiwiii, = (Ax A* x: E (»?^' + ' 

\ i=0 n=0 ^ ^ / 

and the a semi-norm selected for the space of the noise-generating function space 
is: 

\m\B - [Ax At I Y, (r"^' + f^"f 



Claim 32 (New): A method as claimed in claim 16, wherein: 
the a semi-norm selected for the data space is: 

12 


612.42948X00 



and the a semi-norm selected for the space of the noise-generating function space 
is: 


\fi\\B = [Ax At I 


N-i 

E 

n=0 


') 


Claim 33 (New): A method as claimed in claim 17, wherein: 
the a semi-norm selected for the data space is: 



and the a semi-norm selected for the space of the noise-generating function space 
is: 
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Claim 34 (New): A method as claimed in claim 18, wherein: 
the a semi-norm selected for the data space is: 

= (ax At g E + <f] 

\ 7=0 n=0 ^ ^ / 

and the a semi-norm selected for the space of the noise-generating function space 
is: 

(N-l - 
Ax At I 5] 


Claim 35 (New): A method as claimed in claim 19, wherein: 
the a semi-norm selected for the data space is: 
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llulb = ( Ax At £ E ^ + "if) 

\ 1=0 n=0 ^ ^ / 


and the a semi-norm selected for the space of the noise-generating function space 
is: 


m\B = [Ax At I ^ (^-+^ + fr)A 


Claim 36 (New): A method as claimed in claim 20, wherein: 
the a semi-norm selected for the data space is: 

\\,A\i> = [Ax At iZJl-^ + 

\ 1=0 n=0 ^ ^ / 


and the a semi-norm selected for the space of the noise-generating function space 
is: 
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Claim 37 (New ): A method as claimed in claim 21 , wherein: 
the a semi-norm selected for the data space is: 

= [Ax At £ £' ^ + "D' ) 


and the a semi-norm selected for the space of the noise-generating function space 
is: 


Mb = ( Ax At I (r+i + ^r f 


n=0 


Claim 38 (New): A method as claimed in claim 22, wherein: 
the a semi-norm selected for the data space is: 
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and the a semi-norm selected for the space of the noise-generating function space 
is: 


Claim 39 (New ): A method as claimed in claim 23. wherein: 
the a semi-norm selected for the data space is: 

= f Ax A* J: ^ (ur' + u^f) ' 

and the a semi-norm selected for the space of the noise-generating functions space 
is: 
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Claim 40 (New): A method as claimed in claim 24, wherein: 
the a semi-norm selected for the data space is: 


1 



and the a semi-norm selected for the space of the noise-generating function space 
is: 


1 



Claim 41 (New): A method as claimed in claim 25, wherein: 
the a semi-norm selected for the data space is: 
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\\u\\^ = (AxAt + <f] 

\ 1=0 n=0 ^ ^ / 


and the a semi-norm selected for the space of the noise-generating function space 
is: 


Claim 42 (New): A method as claimed in claim 26, wherein: 
the a semi-norm selected for the data space is: 

|l«|lx, = f Ax At g ^ + „r)') 

\ i=0 n=0 ^ ^ / 


and the a semi-norm selected for the space of the noise-generating function space 
Is: 
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Claim 43 (New): A method as claimed in claim 27, wherein: 
the a semi-norm selected for the data space is: 

||»||z> = f Ax At Y/H-^ + ' 

and the a semi-norm selected for the space of the noise-generating function space 
is: 


1 



Claim 44 (New): A method as claimed in claim 1 1 , wherein: 

a distribution of disturbances, in relation to a previously selected reference 
model, of an impedance and of a velocity in the zone of the medium is sought, the 
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correlated noises affecting the data are due to multiple reflections whose kinematics 
and amplitude variations with an offset have been previously estimated, the 
measured data are picked up by seismic surface pickups, the location of the pickups, 
a seismic excitation mode, a recording time and time sampling points being defined, 
and the modelling operator being defined by linearization of the waves equation 
around the model representative of the distribution. 

Claim 45 (New): A method as claimed in claim 44, wherein: 

the cost function quantifying the difference is the square of the semi-norm of 
the difference in the data space. 

Claim 46 (New): A method as claimed in claim 44, wherein: 

adjustment of the model and of the noise-generating functions is obtained by 
a block relaxation method for eliminating unknowns corresponding to each correlated 
noise generating function, the relaxation method being implemented within iterations 
of a conjugate gradient algorithm for calculation of the model. 

Claim 47 (New): A method as claimed in claim 45, wherein: 

adjustment of the model and of the noise-generating functions is obtained by 
means of a block relaxation method for eliminating the unknowns corresponding to 


21 


612.42948X00 


each correlated noise generating function, this the relaxation method being 
implemented within the iterations of a conjugate gradient algorithm for calculation of 
the model. 
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